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Abstract 

The stochastic properties of the total internal energy of a dilute granular gas in the steady 
uniform shear flow state are investigated. A recent theory formulated for fluctuations about the 
homogeneous cooling state is extended by analogy with molecular systems. The theoretical pre- 
dictions are compared with molecular dynamics simulation results. Good agreement is found in 
the limit of weak inelasticity, while systematic and relevant discrepancies are observed when the 
inelasticity increases. The origin of this behavior is discussed. 
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I. INTRODUCTION 



The macroscopic properties of the steady uniform shear flow (USF) of monocomponent 
granular gases have been extensively studied [l-6]. The required quantities to close the 
exact balance equations for the hydrodynamic fields, density, velocity, and internal energy 
or temperature, are the pressure tensor and the cooling rate, since the heat flux vanishes. 
Explicit expressions for all the components of the pressure tensor for arbitrary values of the 
shear rate have been obtained in the literature by using kinetic theory ideas. In addition, 
particle simulation methods, molecular dynamics (MD) and direct simulation Monte Carlo 
(DSMC) techniques, have been employed, and the agreement found between the theoretical 
predictions and the simulation results can be considered as quite satisfactory. On the other 
hand, in all the theories we are aware of, no rheological effects for the cooling rate have been 
reported. Actually, approximate expressions corresponding to zeroth order in the gradients, 
i.e. in the shear rate, are always considered. Although it seems that the cooling rate of 
the steady USF state has never been directly measured in the simulations, the accuracy of 
the predictions for the steady temperature provides strong support for the smallness of the 
corrections to the cooling rate due to the finite shear rate. 

In this work, the focus is put on the stochastic description of a global property of a 
granular gas in the steady USF state, namely the total internal energy. Studying fluctua- 
tions requires using a mesoscopic description of the hydrodynamics of the system, instead of 
the macroscopic description provided by the usual hydrodynamic equations. In the case of 
molecular fluids at equilibrium, the standard mesoscopic theory is the fluctuating hydrody- 
namics proposed by Landau and Lifshitz [zj]. This theory has been extended to systems out 
of equilibrium, mainly in the context of the Navier-Stokes order, although particular states 
beyond it have also been addressed ^illl]. 

Recently, a theory for hydrodynamic fluctuations in an isolated dilute granular gas in 
the homogeneous cooling state (HCS) has been developed by extending the methods of non- 



equilibrium statistical mechanics for molecular gases 15 



17] . The results are expressed as 



Langevin-like equations for the fluctuating density, velocity, and internal energy fields, that 
generalize those by Landau and Lifshitz. The derived equations are not directly applicable to 
other (inhomogeneous) states of a granular gas and, in particular, to the steady USF state. 
Then, what is done in the present work is to formulate, without proof, a generalization of 
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the theory in Ref. [17] based in two main assumptions: 

(i) The stochastic deviations of the hydrodynamic fields from their average values in the 
steady USF state verify the linearization of the general macroscopic hydrodynamic equations 
around that state, supplemented by a fluctuating random force term, i.e. they obey a set of 
coupled linear Langevin equations. 

(ii) The random force terms in the above equations follow from the fluctuating parts of 
the constitutive relations and also from the intrinsic noise induced by the energy dissipation 
in collisions. It is assumed that both noise sources have the same properties as for the HCS, 
with the only difference that the macroscopic quantities associated to the HCS have to be 
replaced by their values in the steady USF state. 

These two assumptions are prompted and stimulated by the results found in molecular 
fluids, where similar properties have been derived by different methods and have proven to 
be useful in many different situations js), [lo|. Since general hydrodynamic equations are only 
available to Navier-Stokes order, the study of fluctuations around the steady USF state has 
to be restricted to that order. Consequently, the results obtained can only be expected to 
hold for small gradients, that for the steady USF state also means small inelasticity, due 
to the coupling between gradients and inelasticity that is characteristic of steady states of 
granular systems. 

The work presented here must be clearly differentiated from other studies in which fluc- 
tuations of driven granular gase s, in contact with some external energy source or thermal 



bath, have been considered 
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141 ] . The latter influences the stochastic properties of the 



granular gas in a nontrivial way and, consequently, it is not clear the relationship between 
fluctuations in granular systems with and without an external source of noise. 

The plan of the paper is as follows. In the next section, the macroscopic properties of the 
steady USF state of a granular gas are shortly reviewed, and the theory of fluctuations around 
that state is formulated. The Navier-Stokes order of approximation is considered, since this is 
the order at which the fluctuating hydrodynamic equations around the homogeneous cooling 
state are known. The application of the theory to the total internal energy of the steady 
USF state is presented in Sec. IIHI The fluctuating energy deviations from the macroscopic 
values obey a linear Langevin equation with a non-white noise. This equation is used in Sec. 
IIVI to compute the energy fluctuations, as well as the two-time correlation function. The 
theoretical results are compared with molecular dynamics simulations of a two-dimensional 
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system of hard disks. Good agreement is found in the quasi-elastic limit, with increasing 
discrepancy as the inelasticity increases, as expected. The final section contains some general 
conclusions and comments. Some of the technical details of the calculations are given in two 
appendices. 



II. FLUCTUATIONS IN THE STEADY USF STATE 



Consider a granular gas composed of smooth inelastic hard spheres (d = 3) or disks 
(d = 2) of mass m and diameter a. The macroscopic balance equations for the number 
density n(r, t), velocity u(r,t), and internal energy density e(r,t), at position r and time t 



arc 
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^ + V-(rm) = 0, (1) 

— + u- Vu + H^V' P = 0, (2) 

at 

— + V • (en) + P : Vu + V • q + (e = 0, (3) 

where P(r,t) is the pressure tensor, q(r,t) is the heat flux, and C( r )^) is the cooling rate, 
the latter following from the energy dissipation in collisions. Equations ([I])-® become a 
closed set of differential equations for the fields n, u, and e, once P, q, and ( are expressed 
as functionals of them by means of the constitutive relations. Then, the equations provide 
a macroscopic or hydrodynamic description of the granular gas, whose domain of validity is 
determined by that of the constitutive relations. 

The steady uniform, or simple, shear flow (USF) state of a granular gas is characterized 
by time-independent and uniform number and internal energy densities, and a velocity field 
given by 



^ij^ji — d^itX^jtyi (4) 



Here, summation over repeated indexes is implicit, Sij is the Kronecker delta symbol, and 
i,j = x, y for d = 2, while i,j = x,y,z for d = 3. Therefore, the velocity profile is linear 
with a constant shear rate a. In addition, the heat flux q vanishes. For this steady flow, 
Eq. ([1]) becomes an identity, Eq. ()2]) implies that the pressure tensor is uniform, and Eq. 
([3]) reduces to 

aP$ + C (0) e(°) = 0. (5) 
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In a hydro dynamic description, Pyx is a given function of a, n^°\ and e^. Consequently, 
Eq. (jSJ) expresses the steady internal energy density as a function of the number density 
and the shear rate. In a compact notation, the hydrodynamic equations following from Eqs. 
©-© plus some constitutive relations can be expressed as 

^ = *,[r, «!{<*>]. (6) 

with [r, t\ {c 7 }] being a nonlinear space functional of the fields {cp} = {n,u,e}. The 
notation indicates that all the space and time dependence on the right hand side of Eq. (J6]) 
is entirely determined from the fields themselves. 

Consider now the fluctuating mesoscopic fields Cp(r, t) whose stochastic average, denoted 
by angular brackets, are the macroscopic fields cp(r,t), 

(Cp(r,t)) = cp(r,t), (7) 

and introduce the deviations 

C'p(r,t) = Cp(r,t)-cf\ (8) 

where Cg' is the value of cp in the steady USF state. Then, it is assumed that the deviations 
C'p obey Langevin-like equations of the form 

r)C (r t) 

-^=Ap[r,t\{cf\C,}] + Mr,t), (9) 
with being a linear term which is related to the hydrodynamic equations by 

Moreover, it is also assumed that the noise terms Tp{r, t) can be identified from the fluctu- 
ating parts of the constitutive relations and the intrinsic fluctuations, both evaluated in the 
HCS. The only change to be made is the replacement of the macroscopic properties of the 
HCS by their values in the steady USF state. 

Equation Q provides a complete formal theoretical scheme to study hydrodynamic fluc- 
tuations in the steady USF state of a granular gas. Nevertheless, a difficulty arises when 
trying to implement it. Although the macroscopic state of the system is known, at least 
approximately, to all orders in the shear rate, to get a similar knowledge of the fluctuations 
and correlations, the most general hydrodynamic equations, i.e. valid for arbitrarily large 



gradient of all the hydrodynamic fields, should be available. Since this is not the case, atten- 
tion will be limited in the following to the Navier-Stokes order. This restricts the expected 
validity of the theory to small values of the shear rate a and, due to the coupling between 
a and a implied by Eq. (J5]) to values of the coefficient of normal restitution close to unity. 

For a dilute gas and keeping only terms up to first order in the gradients of the hydro- 
dynamic fields, the constitutive relations are given by 
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' dui du-j 



OTi a 



(11) 
(12) 
(13) 



q = -ftVT — /iVn, 

where p = 2e/d is the hydrodynamic pressure, T = 2e/nd is the granular temperature, 77 is 
the shear viscosity, k the (thermal) heat conductivity, and /i the diffusive heat conductivity, 
which vanishes in the elastic limit. Moreover, (o is the zeroth order in the gradients cooling 
rate. In the first Sonine approximation, it is 

2C*e 



Co , 
7] = 77*770. 

In the above expressions 770 is the viscosity in the elastic limit 

2 + d 



-r(d/2)7r"V (mT) 



and the dimensionless functions are given by 

2 + d 



4d 



(1"« 2 ) (l 



VV(*-i) 



3a 2 
To 7 



(14) 
(15) 

(16) 



a 2 (a 



16(1 - a)(l - 2a 2 



9 + 24d + (8d - 41)a + 30a 2 - 30a 3 ' 

,n -1 



# , \ _ (3 - 3a + 2d)(l + a) / a 2 
V% ^ ~ Ad V 32 



(17) 

(18) 
(19) 
(20) 



The expressions of the transport coefficients k and fi associated to the heat flux will not be 
needed in the following. 
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III. LANGEVIN EQUATION FOR THE TOTAL INTERNAL ENERGY 



The specific property to be studied in this paper is the fluctuating total internal energy 
of the system. More precisely, the quantity considered is e(t) defined by 

where the integral extends over the volume (for d — 3) or area (for d — 2) V of the system 
and S(r,t) is the fluctuating internal energy density The deviation S'(r,t) = S(r,t) — 
is assumed to obey the corresponding particularization of Eq. ([9]). The equation for 
e(t) follows by integration over the volume (or area) of the system. At this point, the 
boundary conditions play quite a relevant role, as it is also the case even in molecular systems 
at equilibrium. For the purposes of the study being presented here, periodic boundary 
conditions in the Lagrangian frame moving with the local macroscopic velocity field are a 
very convenient choice. First, because they are the boundary conditions consistent with the 
ones used in the numerical simulations that will be reported later on to test the accuracy 
of the theoretical predictions. The second and more fundamental reason is that they lead 
to a decoupling of the equation for e(t) from the equations for the fluctuating number of 
particles and velocity, so that the former obeys a closed simple Langevin-like equation. The 
details of the calculations are given in Appendix while only the result will be reported 
here. The equation obeyed by e reads 

J + Co = C( S ). (22) 

A dimensionless time scale s has been introduced, defined by 

s = — , (23) 

with 

/ e(°) \ 1/2 

v ° = 2 {^m) (24) 

and 

A= (n^V- 1 ) -1 . (25) 

The parameter A is proportional to the mean free path of the gas and v q is a thermal velocity. 
Then, it is easily seen that s is proportional to the average accumulated number of collisions 
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per particle in the time interval between and t. Moreover, C is a dimensionless cooling 
rate defined from by 



AC 



(o) 



v (d + 2)r (d/2) 

The noise term in Eq. (1221) is related with the noise term in the energy equation through 



(26) 



A 



dr J- e (r, t). 



It has two contributions of a rather different physical origin, 

¥ 6 (s) = RJs) + SJs). 



(27) 



(28) 



The properties of the above two noise terms are summarized below. Some details of the 
reasonings leading to them are given in Appendix [B] The contribution R e comes from the 
non-hydrodynamic part of the fluctuating pressure tensor and it has the properties 



(Re(s)) = 0, 



(R e (s)R e (s')) 



4Co 
Ndrj 



G(\s-s'\). 



(29) 



(30) 



Here N = n^V is the total number of particles, r\ is the dimensionless shear viscosity 
defined as 



'I 



Tj 



(0) 



n^mXvo 



and the function G(s) has the form 



with 



G(s) 



Ka) 



1 + 02(a) ,T 



4 



A 4 = ((a) + 



41(a) 



(2d + 3-3a)(l + a)7r— 



1 + 



1 + a 2 (a 
23a2(a) 
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(31) 



(32) 



(33) 



2\/2d(d + 2)r (d/2) 

The second noise contribution on the right hand side of Eq. fl28l) arises directly from the 
fluctuations of the energy dissipation in phase space, and vanishes in the elastic limit a — > 1. 
It is 

(Se(s)) = 0, (34) 

and 

(S £ (s)S e (s')} = ^( a 33 (a)5(s - s'), (35) 
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where 



6(a) 



. . d+1 d + 2 . . . 
033(a) = + —j-j- 02(a) + 6(a), 

2 + d - 6d 2 - (10 - 15d + 2d 2 )a - 2(2 + 7d)a 2 + 2(10 - d)a 2 



Qd(2d + 1) - 2d(ll - 2d)a + 12da 2 - 12da 3 
Moreover the two noise terms are uncorrelated, 

(R e (s)S e (s')) = 0, 



(36) 
(37) 

(38) 



for all s and s'. Also, as usual in the Langevin description, the fluctuating noise term at a 
given time is not correlated with the energy variable at previous times, 



(MsHs')) = 0, 



(39) 



for s > s'. 



IV. ENERGY FLUCTUATIONS AND TIME-CORRELATION FUNCTION 



The general solution of Eq. ([22]) reads 

e( s ) = e~^°e(0) + [ ds'e- {s - s,) ^ e (s') 
Jo 

Taking into account Eq. (|39|) . it is found that 



(40) 



(e 2 ( s )) = e - 2 ^(e 2 (0))+ / ds' ds"e-V s - s '- s "K°(® e (s')$ e (s")). (41) 

Jo Jo 

For times s ^> (q 1 , the first term on the right hand side of the above equation can be 
neglected, and using Eqs. ( |30l) . ( 1351) . and ( )38l) . one gets, 

€0 



N 



( e 2 (s)) = ^ / ds > / rfs " e -(^-'-")Co 



x 



— G(| S '- S "|) + «33(a)o( S '- S A 
77a 



where the function G(s) is defined in Eq. fl32|) . An easy calculation leads to 



drG{T)e- T ^ -e~ 2sC ° / rfrG(r)e r ^ 
IV Jo 



(42) 



(43) 



Considering again the limit s 3> Co an d keeping in mind that G(r) decays faster than 
exp(— rCo)j the steady value 



(e 2 ) = - 



2a 33 (a) + — / drG(r)e- rC « 
Vd Jo 



(44) 



is obtained. Finally, evaluation of the integral on the right hand side yields 

1 + 2a 2 (a) 



(e 2 > = - 



2a 3 3(a) + 



477|I(a)|dJ ' (45) 
A term quadratic in a2(a)has been neglected by consistency with the approximation in 



which Eq. ffl8|) is derived 
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23]. 



To check this theoretical prediction, Molecular Dynamics (MD) simulations of a system 
of inelastic hard disks in a square box of side L have been performed, using Lees-Edwards 
boundary conditions 24| . These are periodic boundary conditions in the Lagrangian frame 
moving with the local velocity of the fluid, so that in that frame the system is macroscopically 



lomogeneous 



The simulations have been carried out using an event driven algorithm 



261 ] . The number of particles in the system was A" = 2000, corresponding to a number 



density = 0.02cr -2 . In all the simulation results to be reported in the following, it 
was checked that the system reached, after a transient period, a steady state with the 
macroscopic profiles corresponding to the USF. The steady values shown in the following 
have been obtained by averaging over a number of trajectories, typically 100, and also on 
time, about 150 collisions per particle. 

In Fig. [1], the results obtained for the second moment of the normalized internal energy 
fluctuations, (e 2 ) are plotted as a function of the coefficient of normal restitution a. The 
data shown have been obtained in a system with a = 6, 32 x 10~ 3 (T(0)/m) 1 ^ 2 cr -1 , where 
T(0) is the initial temperature, but equivalent simulation results have been obtained with 
other values of the shear rate, in agreement with Eq. (145!) that does not depend on a. 
This function is also plotted in the figure (solid line). In the small inelasticity region, say 
a > 0.95, a fairly good agreement between theory and simulation is observed, but when the 
coefficient of restitution is further decreased, significant quantitative discrepancies show up. 
This behavior was expected, since the Navier-Stokes approximation being used here implies 
small inelasticity as discussed above. Let us mention that the shape of the probability 
distribution for the internal energy fluctuations has also been investigated, and that it is 
always very well fitted by a Gaussian, within the statistical uncertainties. 

A trivial way of incorporating some rheological effects into the theory is to substitute in 
Eq. (1451) the Navier-Stokes shear viscosity by a generalized shear viscosity, valid to all orders 
in the gradients, and defined by 

(o) 



(0) Pxy . 

Vg = — — ■ ( 46 ) 

(Jj 

10 



2.5 



2 

<£ > 



1.5 



A 




FIG. 1: (Color online) Dimensionless second moment of the total energy fluctuations (e 2 ) of a 
granular gas of hard disks in the steady USF state, as a function of the coefficient of normal 
restitution a. The (red) symbols are from MD simulations of a system of hard disks of density 
n(°) = 0.02cr~ 2 , the solid line is the Navier-Stokes theoretical prediction given by Eq. (|45p . and the 
dashed line has been obtained by using the rheological expression for the viscosity instead of the 
Navier-Stokes one. 



Several theoretical expressions for the 



have been obtained in the literature 



0£ 



generalized viscosity T]q of a dilute granular gas 



271 ] . They are practically indistinguishable in the 



range of values of a considered in Fig. [T] When any of these expressions is used into Eq. 

after reducing it accordingly to Eq. (15T|) . the result indicated in Fig. [T]by means of 
the dashed line is obtained. Although there is a significant improvement as compared with 
the Navier-Stokes approximation result, it is clear that this effect is not enough to explain 
the observed increase of the energy fluctuations as the coefficient of restitution decreases. 
The time correlation function Ce(£, if) for the total energy is defined as 

C E {t,t') = JdrJ dr'(S(r,t)£(r',t')) - V 2 e m = V 2 e m (e{t)e{t')), (47) 
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for t < f < 0. Use of Eq. ([22]) gives 

( e (t)e(t')) = e-^o(e 2 ), (48) 

with 

s-s'= Vo{t - f \ (49) 
A 



In the derivation of the above results, Eq. (1391) has been employed. Therefore, the theory 
of fluctuations developed here leads to the rather strong prediction that the decay of total 
internal energy fluctuations of a dilute granular gas is governed by the same quantity in the 



steady USF state as in the homogeneous cooling state [15], being exponential in both cases. 

Results obtained from MD simulations for the decay of the energy time correlation func- 
tion are shown in Fig. [2] using a logarithmic scale. The plotted quantity is C^(s, s') = 
C 'E(t,t') I 'C 'e(1? ,t') . An exponential decay is clearly identified, since the deviations observed 
for large times are due to lack of statistics. 

By fitting the decay of the time correlation function Ce to an exponential, its decay rate 
has been estimated as a function of a. The theoretical prediction is that it agrees with the 
zero order in the gradients cooling rate ( (see Eq. ( 148]) ). The comparison is carried out in 
Fig. [3] Again, a good agrement is found for small inelasticity and an evident systematic 
discrepancy shows up as the coefficient of restitution decreases. It is worth to remark that 
in this case the theoretical prediction can not be improved by incorporating in a trivial way 
rheological effects, as done above for the energy fluctuations. It has already been indicated 
that there are strong evidence that rheological corrections to the cooling rate, if any, are not 
relevant. The consequence seems to be that some new physical effects must be incorporated 
into the theory to extend it to the not small inelasticity region. 

V. DISCUSSION 

In this paper, a recent theory of fluctuations and correlations for a dilute granular gas 
in the homogeneous cooling state 17| has been extended in a nontrivial way to a system 
in the steady uniform shear flow state. The extension has been carried out by similarity 
to what happens in molecular systems. From the theory, predictions have been derived for 
the second moment of the internal energy fluctuations and also for its two-time correlation 
function. They have been compared with molecular dynamics simulation results, and good 
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FIG. 2: (Color online) Normalized dimensionless internal energy time correlation function C* E as 
a function of the time interval s — s' for several values of the coefficient of normal restitution, as 
indicated in the inset. This dimensionless time interval is proportional to the accumulated number 
of collisions per particle. The density is the same as in Fig. [TJ 

agreement has been found in the small inelasticity region, while discrepancies become evident 
and systematic as the inelasticity increases. This was expected, as the theory is formulated 
at the level of the Navier-Stokes order of approximation, and limitation to small gradients 
also implies limitation to small inelasticity, because of the coupling between both parameters 
in the steady USF state. Nevertheless, there is a conceptually relevant difference between 
fluctuations and time correlations. While in the former, the role played by rheological effects 
can be easily foreseen, for instance through the viscosity, in the latter non-Newtonian effects 
seem to be negligible for the cooling rate, that is the quantity governing the rate of decay 
to Navier-Stokes order. Therefore, some new physical effects are needed to explain the 
failure of the theory to describe the decay of the energy time correlation function for high 
and moderate inelasticity. In particular, this requires to consider the transport coefficients 
characterizing the response of the system to small perturbations of the steady shear state 



13 




FIG. 3: (Color online) Dimensionless decay rate of the internal energy time correlation function in 
the steady USF state. The symbols are numerical results obtained by fitting the results in Fig. [2j 
The solid line is the theoretical prediction, i.e. the zeroth order cooling rate given by Eq. (|26p . 
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29]. 



It is worth to stress the two main differences of the Langevin equation for the total 
internal energy of the steady USF state as compared with similar equations for equilibrium 
molecular states. Firstly, the noise source term coming from the fluctuating part of the 
pressure tensor is not delta correlated in time, i.e. the noise is not white. Secondly, there 
is another noise term that is intrinsic to the energy dissipation and vanishes in the elastic 
limit. Both terms are relevant, in the sense of leading to non-negligible contributions to the 
energy fluctuations and time correlations. This is an indication that extending the theory 
of equilibrium fluctuations to non-equilibrium situations can be far from trivial. 
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Appendix A: Linearization of the Navier-Stokes equations around the steady USF 
state 

Define the fluctuating fields 

Af (0) (r,t) =Af(r,t) -n {0 \ (Al) 

W ,(0) (r,t) =U(r,t) -u (0) , (A2) 

S'(°\r,t) =S(r,t)-e {0 \ (A3) 

where Af(r,t), U(r,t), and S(r,t) are the fluctuating number density, velocity field, and 
internal energy density, respectively. The aim of this Appendix is to derive the equation 
for e(t), defined in Eq. (I2TI) . Accordingly with the theory formulated in this paper, this 
requires to linearize the hydrodynamic Navier-Stokes equation for the energy by applying 
Eq. ([nD followed by a space integration over the volume or surface of the system. In the 
process, a relevant role will be played by the assumed periodic boundary conditions for the 
fluctuations of the hydrodynamic fields C'l , (r,t). Each of the terms in the Navier-Stokes 
equation for the energy will be now considered separately. Application of the linearization 
indicated in Eq. (jUJ) to the convective term gives 

V • (en) — > e (o) v . W '(») +u (»).V£'(°). (A4) 

Upon writing the above equation the form of the macroscopic velocity field of the USF has 
been taken into account. Moreover, it is 

J drV-U m (r,t) = (A5) 

and 

J dr n (0) ■ V£ /(0) (r, t) = j dr <«/ J^' (0) (r, t) = 0, (A6) 

because of the periodic boundary conditions. Consequently, there is no contribution to the 
equation for e{t) due to the convective term of the hydrodynamic equation. Actually, this is 
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a general result, that is not restricted to the Navier-Stokes order. Next, consider the term 
involving the pressure tensor in Eq. ([3]). It leads to 

P:V«^ P (0) : VW ,(0) + aV'$\ (A7) 

where Vyx is the deviation of yx component of the fluctuating pressure tensor from its 
macroscopic value in the steady USF state. Due to the homogeneity of P^ ^ it follows by 
using Eq. ( 1A5I) that the first term on the right hand side vanishes when integrated over 
the system. To analyze the term proportional to the fluctuation of the pressure tensor, the 
constitutive Navier-Stokes relations will be taken into account. Then, 

+ Mw) <v,(r ' () ' (A8) 

Use of Eq. <^2) yields 

f^|-( W ,-,<), (A9) 

because rj oc (e/n) 1 ^ 2 . Then, 

/*•/ *' (^iT^''^ drM ' te ^^ = °- < A10 > 

again because of the periodic boundary conditions and the conservation of the number of 
particles. For the same reason it is 

jdrj dr> ( ^P| ) (0) ^ (0) K^)--^ (0) / dr [V^(r, t) + V x U^(r, t)] = 0. 

(All) 

The last contribution from the fluctuating pressure tensor involves 
so that, 

*' {WwT )£ ' i ° )<r ' j) = -^'"^ < ai3) 

This is the only non- vanishing contribution arising from the linearization of the viscous term. 
The heat flux term gives 

V • q -> V • Q, (A14) 
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where Q is the fluctuating heat flux. Although its explicit expression to Navier-Stokes order 
can be easily obtained, it is not needed here, since the periodic boundary conditions directly 
imply 

drV-Q = 0. (A15) 



This result holds beyond the Navier-Stokes order and it is valid to all orders in the gradients 
of the hydrodynamic fields. 

The last term to be linearized is the one describing the energy dissipation in collisions, 



Coe->Cr^ (0) + eo4 (0) - 



(A16) 



The fluctuating cooling rate Z'(°> is 



^"{WaT**™ (Al7) 



(A18) 



The only term leading to a non-zero contribution after integration over r is 

where it has been used that £o oc (ne) 1 ^ 2 . Putting together all the above results, it follows 
that e(t) obeys the evolution equation 



at + e 



e(t) = * e (t), 



where 



and 



drT e (r,t) 



(A19) 



(A20) 



„(0)„2 o/(0) 
O (0) _ t ' a , - a(o) 

£ ~ 2e(°) + 2 

Upon writing the last equality, it has been used that to Navier-Stokes order Eq. (j^J) reads 



2e( o)+^ = Cr. (A21) 



y>) _ c (o) e (o) 



0. 



(A22) 



Equation (1221) follows directly from Eq. (1A19I) by making the change of time scale given in 
Eq. (1231). 
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Appendix B: The noise term in the Langevin-like equation for the internal energy 



Let R(r, s) be the non-hydrodynamic part of the fluctuating pressure tensor. From the 
energy balance equation, Eq. ([3]), it is seen that 

Re(s) = — ^— / dr Ry X (r,s). (Bl) 



Ve(°)v 

In ref. |l7|, a noise term R^ MG (k, s) was considered (see Eqs. (53) and (54) there) for the 
HCS. Its definition is related to Rij(r,s) by 

^ e A ->BMG 



drR yx (r,s) = —^—R^(k = 0,s). (B2) 
The factor on the right hand side arises because of the scaling of the length and velocity in 



ref 



17| . Then, assuming that the results for the HCS can be mapped into relations for the 



steady USF state as discussed in the main text, Eqs. ( )29l) and ( !30|) here are easily derived 
from Eqs. (53) and (54) in fl7|. 

The second term on the right hand side of Eq. ( 1281) is a direct consequence of the 
inelasticity of collisions. Denoting by S e (r,t) the intrinsic noise term appearing in the 
equation for the fluctuating energy £'^(r,t), it is 

Se(s) = -^r- f drS e (r,s). (B3) 



VeWvo 

This intrinsic noise term has also been studied in ref. [17[ for the HCS. It is easily verified 
that 

J dr S e (r, s) = \ d - l e {0) v Q Sf MG {k = 0, s), (B4) 

where Sf MG (k, s) is the intrinsic noise term defined in [l?]]. Properties ( 13~4"|) and (I3"5"|) in this 
paper are then equivalent to properties (29) and (D.4) in that reference. Of course, this 
reasoning requires to assume that the properties of the intrinsic noise in the HCS are the 
same as in the steady USF state, the only difference being in the macroscopic fields. 
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